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Abstract 

A  field  deployable  Tunable  Diode  Laser  Atmospheric  Spectroscopy  (TDLAS) 
device  has  been  used  to  investigate  the  atmospheric  absorption  of  Diode  Pumped  Alkali 
Laser  (DPAL)  wavelengths  at  a  distance  of  2,000  meters  in  a  maritime  environment.  The 
spectral  regions  surrounding  the  cesium  (Cs)  and  rubidium  (Rb)  DPAL  emission  lines, 
890  to  910  nanometers  and  790  to  800  nanometers  respectively,  were  examined  in  order 
to  determine  the  effects  of  absorption  by  the  (0  0  0)— >  (0  0  3)  and  (0  0  0) — >(  1  0  2)  water 
vapor  vibrational  bands. 

Spectral  data  for  the  Cs  spectral  region  were  analyzed  to  determine  temperature 
estimates  of  between  6.60  to  30.02°  C  with  statistical  errors  of  between  0.1  and  1.14°  C 
and  concentration  estimates  of  between  2.42  to  3.5 1  x  10  cm"  with  error  bounds  0.01 1 
and  0.001  x  10  cm"  over  the  9  day  span  of  data  collection.  Temperature  and  pressure 
estimates  were  within  6%  and  25%,  respectively,  of  those  reported  by  the  Davis  Vantage 
Pro2  portable  weather  station  used  to  collect  on  site  meteorological  data  for  the 
experiment. 

The  observed  spectra  were  compared  to  data  previously  collected  at  150  and 
1,000  meters  and  the  differences  were  investigated.  The  strength  of  the  absorption 
features  at  the  150  and  2,000  meters  scaled  primarily  with  path  length  as  the  water  vapor 
concentration  only  varied  slightly  between  the  two  experiments.  The  strengths  at  1,000 
and  2,000  meters  scaled  in  part  with  path  length  and  in  part  with  concentration,  with  the 
combination  of  these  effects  only  differing  from  the  actual  intensity  ratio  by  15%. 


IV 
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The  rotational  dependence  of  broadening  rates  was  investigated.  By  analyzing  the 
widths  of  the  [1  1  0] — >[2  1  1]  and  [6  2  5]— >  [7  2  6]  rotational  transitions  in  the  (0  0 
0)— >(0  0  3)  vibrational  band  of  water  vapor,  it  was  found  that  the  observed  linewidths 
were  approximately  25%  higher  than  those  reported  in  the  High  resolution 
TRANsmission  molecular  absorption  (HITRAN)  2008  database  and  that  the  scaling  ratio 
of  the  observed  rates  was  5%  below  the  ratio  of  the  reported  values. 
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INVESTIGATION  OF  THE  ATMOSPHERIC  PROPAGATION  OF  ALKALI  LASERS 


IN  A  MARITIME  ENVIRONMENT  USING  TUNABLE  DIODE  LASER 
ATMOSPHERIC  SPECTROSCOPY 

I.  Introduction 

Significant  effort  into  the  development  of  Diode  Pumped  Alkali  Lasers  (DPAL’s) 
began  in  2003  with  Krupke’s  demonstration  of  the  first  3  level  rubidium  (Rb)  laser  [1], 
After  this  first  demonstration,  which  used  a  titanium  sapphire  laser  as  the  pump  source, 
work  began  on  semiconductor  laser  pumped  systems,  resulting  in  the  first  demonstration 
of  a  true  DPAL  by  Page  in  2006  [2].  Since  then,  the  pace  of  research  efforts  regarding 
alkali  lasers  have  increased  significantly,  resulting  in  similar  cesium  (Cs)  and  potassium 
(K)  laser  demonstrations  in  2003  and  2007  respectively  [3,  4].  Subsequent  efforts  in 
improving  all  three  of  these  systems  have  resulted  in  a  kW-class  Cs  system  with  a 
conversion  efficiency  of  48%  [5]. 

DPAL’s  have  attracted  the  attention  of  the  laser  community  for  several  reasons. 
Due  to  the  small  spectral  spacing  between  the  upper  two  levels  of  any  3  level  DPAL, 
about  15  nm  for  Rb  and  5  nm  for  Cs,  these  systems  have  quantum  efficiencies  of  between 
95  and  99%.  This  allows  mixing  between  the  upper  two  states  via  a  buffer  gas  and 
reduces  the  total  amount  of  waste  heat  produced  in  the  gain  medium.  In  addition,  the  fact 
that  the  gain  medium  is  a  vapor  makes  overall  thermal  management  of  the  gain  medium 
easier  than  in  a  solid  state  laser  system. 
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When  used  in  conjunction  with  a  suitable  diode  laser,  a  DPAL  can  be  used  as  an 
electrically  driven  system  with  no  need  for  hazardous  chemicals  or  complicated  pump 
systems.  In  addition,  due  to  their  exceptionally  high  gain  coefficients,  orders  of 
magnitude  higher  than  that  of  other  gas  phase  gain  media,  it  is  possible  for  a  significant 
amount  of  power  to  be  extracted  from  a  small  gain  volume.  These  factors,  combined 
with  an  exceptional  beam  quality  due  to  the  gaseous  nature  of  the  gain  medium,  result  in 
a  system  that  has  come  under  investigation  for  use  in  high  energy  laser  applications. 

As  part  of  this  investigation,  it  is  important  to  characterize  atmospheric  absorption 
in  the  near-IR  region.  For  example,  it  is  known  that  several  water  vapor  absorption  lines 
surround  both  the  Cs  and  Rb  DPAL  emission  lines.  The  spectral  region  around 
potassium  (K),  however,  contains  several  strong  oxygen  (CL)  absorption  features.  These 
absorption  features  will  affect  a  propagating  laser  beam  based  on  their  location  with 
relation  to  the  emission  line  as  well  as  the  strength  of  the  absorption  feature. 

Due  to  the  20  nm  tunable  range  of  the  diode  heads  used,  it  was  possible  to  scan 
over  the  majority  of  the  R-branch  of  the  (0  0  0)— >(0  0  3)  FLO  vibrational  band.  By 
observing  the  width  of  two  of  the  rotational  lines  in  this  branch,  it  was  possible  to 
investigate  the  rotational  dependence  of  broadening  rates.  These  values  were  then 
compared  to  broadened  half  widths  reported  in  the  HITRAN  2008  database. 

Previous  TDLAS  investigations  at  DPAL  wavelengths  have  been  conducted  at 

path  lengths  of  150  m  and  1  km.  By  collecting  absorption  spectra  for  a  2  km  path,  it  is 

possible  to  more  accurately  predict  how  a  high  energy  laser  (HEL)  system  using  a  DPAL 

would  perform  over  longer  distances.  This  comparison  also  allowed  for  an  investigation 

of  the  dependence  of  absorbance  on  path  length  and  concentration.  Doubling  the  path 
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length  allows  for  a  better  understanding  of  the  difficulties  associated  with  performing 
future  TDLAS  investigations  at  increased  path  lengths.  In  addition,  this  test  was 
performed  in  a  maritime  environment,  allowing  for  an  investigation  of  how  a  higher 
concentration  of  water  vapor  affects  the  propagation  of  a  DPAL.  Characterizing  the 
extent  to  which  these  absorption  features  degrade  the  performance  of  an  alkali  laser  at 
longer  path  lengths  and  in  different  physical  locations  is  an  important  step  in  determining 
the  suitability  of  a  system  for  various  HEL  applications. 
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II.  Background 


TDLAS  Background 

While  conventional  spectroscopy  uses  a  broadband  source  and  a  resolving  device 
such  as  spectrometer  or  interferometer  to  interrogate  a  sample,  laser  spectroscopy  uses  a 
spectrally  narrow  tunable  laser.  Through  the  use  of  a  sufficiently  narrow  laser  source, 
such  as  a  dye  laser  or  narrowed  diode  laser,  it  is  possible  to  measure  absorption  features 
with  a  resolution  limited  by  the  linewidth  of  the  atomic  or  molecular  transition  of  interest. 
The  introduction  of  tunable  diode  lasers  makes  laser  spectroscopy  much  more  convenient 
by  cutting  out  the  more  traditional  complex  dye  laser  systems  and  replacing  them  with 
electrically  driven  systems.  While  they  typically  have  lower  output  powers  than  dye 
lasers,  tunable  diode  lasers  offer  a  linewidth  of  less  than  1  MHz  and  can  typically  be 
tuned  over  a  range  several  nanometers.  These  performance  characteristics,  in  addition  to 
their  reduced  complexity  and  size  compared  to  dye  laser  systems  make  TDLAS  an 
excellent  tool  for  investigating  molecular  absorption  in  the  atmosphere  with  resolution 
limited  by  pressure  broadening. 

H20  Spectroscopy 

Because  of  its  bent  triatomic  structure,  a  water  molecule  has  three  vibrational 
modes.  The  first,  Vi,  is  referred  to  as  a  symmetric  stretch.  In  this  mode,  both  hydrogen 
atoms  oscillate  back  and  forth  with  respect  to  the  oxygen  atom  symmetrically.  The 
second,  v2,  is  referred  to  as  a  bend.  In  this  mode,  the  molecule  straightens  out  and  flexes 
again.  The  third  mode,  V3,  is  the  asymmetric  stretch.  In  this  mode,  one  hydrogen  moves 
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toward  the  oxygen  molecule  while  the  other  moves  away.  These  three  modes  make  up 
the  vibrational  quantum  number  for  H2O  and  are  typically  represented  as  (vi,  V2,  V3).  [6] 
A  diagram  of  each  of  these  modes  is  shown  in  Figure  1  below. 


Figure  1:  Three  vibrational  modes  of  a  water  molecule,  the  symmetric 
stretch  (left),  the  bend  (center),  and  the  asymmetric  stretch  (right). 

The  rotational  states  within  a  vibrational  band  are  represented  in  a  similar  way. 

Because  water  is  an  asymmetric  rotor,  each  rotational  transition  must  be  represented  by 

three  rotational  quantum  numbers,  J,  Ka,  and  Kc.  The  primary  rotational  quantum 

number,  J,  represents  the  total  amount  of  rotational  energy  contained  by  the  molecule, 

with  Ka  and  Kc  representing  the  two  symmetric  limits  of  the  possible  orientations  of  the 

hydrogen  atoms  with  respect  to  the  oxygen  atom  which  taken  together  define  how 

asymmetric  the  molecule  is.  These  three  quantum  numbers  are  typically  listed  as  [J,  Ka, 

Kc],  similar  to  the  vibrational  quantum  numbers  [6].  Using  these  two  representations,  it 

is  possible  to  characterize  a  spectrum  using  the  vibrational  band  and  the  rotational 

transition  within  that  band. 

Molecular  Absorption 

The  primary  principle  at  work  in  TDLAS  is  Beer’s  Law.  Beer’s  Law  relates  the 
ratio  of  the  intensity  transmitted  through  a  sample,  7t,  and  the  incident  intensity,  Iq,  with 
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the  absorption  cross  section,  a,  species  concentration,  N,  and  path  length,  /  as  a  function 
of  the  frequency  in  wavenumbers,  v,  by  the  relation 

it(v)  =  i0(y)e-*w-Nl  (l) 

with  the  absorption  cross  section  given  by  the  relation 

=  f^'WCv)  (2) 

where  g2  and  gj  and  the  degeneracies  of  the  two  states  of  interest,  X  is  the  wavelength  of 
the  transition  at  the  line  center,  v  is  frequency  of  interest  expressed  in  wavenumbers  or 
inverse  centimeters.  A21  is  the  spontaneous  emission  coefficient,  and  g(v)  is  the  lineshape 
function,  taking  the  form  of  a  Gaussian,  Lorentzian,  or  Voigt  profile.  At  times,  it  is  more 
convenient  to  combine  the  absorption  cross  section  with  the  species  concentration  in  one 
variable,  the  absorption  coefficient,  a(v)  which  can  be  used  to  rewrite  Beer’s  Law  as 

It(v)  =  I0(y)e~am-1.  (3) 


Lineshapes 

The  lineshape  function,  g  (v) ,  can  have  one  of  three  forms.  The  first,  the 
Lorentzian  profile,  takes  the  form 


9l(v) 


AvL/(2n) 

(^)2+(*-*o)2' 


(4) 


where  v0  is  the  frequency  at  line  center  expressed  in  terms  of  wavenumbers  and  AvL  is 
the  FWHM  of  the  a  Lorentzian  profile  in  wavenumbers,  given  by 


AVl  2 PtoTT^ik 


nfikT 


(5) 
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where  pt0rr *s  the  pressure  in  torr,  olk  is  the  collisional  cross  section  for  species  of 
interest,  p  is  the  reduced  mass  of  the  collision  pair,  k  is  the  Boltzmann  constant,  and  T  is 
temperature  in  Kelvin.  The  lineshape  takes  this  form  in  the  case  of  homogenous  pressure 
broadening,  meaning  that  as  you  increase  the  pressure,  the  line  will  broaden  regardless  of 
its  spectral  relation,  typically  at  a  rate  of  around  10-20  MHz/torr.  In  the  case  of  open 
path  TDLAS,  the  pressure  is  around  1  atm,  which  equates  to  760  torr.  This  results  in  a 
typical  pressure  broadened  linewidth  of  -10  GHz,  much  larger  than  the  instrument 
limited  resolution  dictated  by  the  laser  source  of  ~300MHz. 

The  second  form  the  lineshape  function  could  take  is  that  of  a  Gaussian,  the  form 
of  which  is  given  as 


In  (2) 


(6) 


where  AvD  is  the  FWHM  of  a  Doppler  broadened  lineshape,  given  as 

Avd  =  7.2x10  ~7v0J^.  (7) 

where  T  is  temperature  in  Kelvin  and  M  is  mass  in  atomic  mass  units.  The  Gaussian 
lineshape  is  found  in  the  case  of  inhomogeneous  broadening  caused  by  an  increase  in 
temperature  [7].  For  typical  values  of  temperature  experienced  over  the  course  of  this 
experiment  of  around  290  K,  the  FWHM  is  ~1  GHz,  again  larger  than  the  instrument 
limited  resolution.  However,  this  is  an  order  of  magnitude  below  the  pressure  broadened 
FWHM,  meaning  that  the  primary  factor  governing  the  overall  broadening  rate  of 
spectral  lines  is  pressure. 
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While  it  is  an  order  of  magnitude  higher  than  the  Doppler  broadened  width,  the 
pressure  broadened  FWHM  is  not  sufficiently  large  to  entirely  dominate  the  overall 
FWHM.  Because  of  this,  it  is  important  to  also  consider  the  temperature  broadening  rate 
when  analyzing  the  overall  lineshapes  in  the  data  collected  over  the  course  of  this 
experiment.  This  is  done  through  the  use  of  a  Voigt  profile,  a  convolution  of  the  relevant 
Lorentzian  and  Gaussian  lineshapes.  The  HITRAN  1996  documentation  provides  a 
method  for  using  parameters  found  in  the  HITRAN  database  to  determine  the  half  width 
of  an  absorption  feature  through  the  relation 

y(.P>  T)  =  (pf-)  C Yair{prefJref)iV  ~  Vs )  +  Yself{PrefJref)Ps)-  (8) 

where  T  and  Tref  are  the  actual  and  reference  temperatures  in  Kelvin,  p,  ps,  and  pref  are  the 
actual  pressure,  partial  pressure  of  the  molecule  of  interest,  and  reference  pressure  in 
atmospheres,  ywr  and  yseif  arc  the  air-broadened  and  self-broadened  half  widths  found  in 
the  HITRAN  database,  and  n  is  the  exponent  dictating  the  temperature  dependence  of  the 
transition,  also  found  in  the  HITRAN  database.  y(p,  T)  is  the  actual  half  width,  corrected 
from  the  reported  half  widths  to  account  for  temperature  and  pressure  [8]. 

The  partial  pressure  of  the  molecule  of  interest,  in  this  case  water  vapor,  is  found 
by  multiplying  the  vapor  pressure  of  water  in  the  atmosphere  by  the  relative  humidity  in 
the  atmosphere.  The  vapor  pressure  was  found  using  the  relation 

Pvap  =  iog10(Pst )  +  8.1328  X  10-3  (io-3-49149^-1)  _  i)  +  io'7-90298^-!)  + 
5.02808  log10  ~  13816  x  10~7  (9) 
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where  pvap  is  the  vapor  pressure  in  hPa,  Tst  is  the  vaporization  temperature  of  water,  T  is 
the  actual  temperature,  and  pst  is  the  1  atmosphere  in  hPa  [9], 

Water  vapor  exhibits  over  250  absorption  lines  over  the  890  nm  to  910  nm  range 
surrounding  the  Cs  DPAL  emission  line.  These  lines  are  described  by  independent  Voigt 
lineshapes,  which  are  superimposed  onto  another  on  a  frequency  axis  to  form  a  complete 
spectrum.  With  information  about  precise  transition  frequencies  and  absorption 
coefficients  for  each  of  these  independent  lines  and  knowledge  of  the  temperature, 
pressure,  and  concentration  during  an  experiment,  it  is  possible  to  use  these  relations  to 
generate  a  spectral  model  with  which  to  compare  experimental  results. 

The  reverse  is  also  true.  Using  the  Line  By  Line  Radiative  Transfer  Model 
(LBLRTM)  and  the  parameters  found  in  the  HITRAN  database,  it  was  possible  to 
simultaneously  fit  a  Voigt  profile  to  each  known  spectral  line  in  a  collected  data  set  using 
temperature,  pressure,  and  concentration  as  variable  parameters  in  order  to  obtain 
estimates  for  those  variable  parameters  listed.  In  the  absorption  coefficient,  the  species 
concentration  is  found  multiplied  by  the  absorption  cross  section.  This  means  that  an 
increase  in  concentration  will  have  the  effect  of  increasing  the  strength  of  every 
absorption  line  in  the  spectrum  by  the  same  ratio.  Atmospheric  pressure  is  found  within 
the  Lorentzian  contribution  to  the  Voigt  profile.  As  pressure  increases,  the  lineshape  will 
broaden,  meaning  that  the  majority  of  pressure  information  is  contained  in  the  wings  of 
an  absorption  feature.  Unlike  concentration  and  pressure,  temperature  is  determined  not 
by  Beer’s  Law  parameters,  but  rather  by  the  Boltzmann  intensity  distribution  of  the 
rotational  energy  levels  in  the  absorption  spectrum.  As  the  temperature  changes,  the 
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rotational  transition  with  the  highest  intensity  will  change.  The  peak  rotational  transition 


is  found  using  the  relation 


where  k  is  the  Boltzmann  constant,  T  is  the  temperature,  and  Bv  is  the  rotational  energy 
constant  for  a  given  vibrational  transition.  With  a  rotational  constant  of  435.346  GHz  for 
water  and  an  atmospheric  temperature  of  290  K,  the  calculated  maximum  rotational 
quantum  number  is  2.134.  This  compares  favorably  with  the  observed  peak  of  ~2.  All 
other  parameters  which  determine  lineshape,  such  as  species  mass,  line  center  location, 
and  stimulated  absorption  coefficients  are  fixed  values  for  a  given  species  and  will  not 
change  over  the  course  of  an  experiment. 

Molecular  Scattering 

In  addition  to  molecular  absorption,  Rayleigh  scattering  also  affects  the  way  in 
which  light  propagates  through  the  atmosphere.  This  is  the  scattering  of  incident  light 
due  to  particles  much  smaller  than  the  wavelength  of  that  light,  typically  atoms  or 
molecules.  Rayleigh  scattering  has  a  A"4  dependence,  meaning  that  shorter  wavelengths 
are  scattered  more  than  longer  wavelengths.  Because  of  this  dependence,  scattering 
introduces  a  baseline  curvature  of  around  9%  from  890  nm  to  910  nm.  This  scattering  is 
independent  of  absorption,  meaning  that  it  is  possible  for  scattering  to  obscure  absorption 
features  if  baseline  noise  is  not  the  dominant  effect.  The  scattering  coefficient  is 
typically  on  the  order  of  10"  cm"  ,  meaning  that  scattering  effects  alone  can  conceal 

2  5 

absorption  features  lower  than  10  "  in  absorbance  at  a  distance  of  10  cm.  This  limit  is  an 
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order  of  magnitude  lower  than  the  observed  noise  in  the  baseline  of  the  data  collected  at 


2,000  m. 


Turbulence 

Atmospheric  turbulence,  when  significant,  can  also  degrade  the  performance  of  a 
laser  at  long  propagation  distances.  The  small  changes  in  temperature  along  the  path  of 
the  beam  induce  changes  in  the  index  of  refraction,  which  bend  and  distort  the  beam  as  it 
travels.  The  resulting  expanded  effective  diameter  of  the  distorted  beam  a  distance  R 
away  from  the  laser  aperture  is  given  by  the  relation 

V2 1R 


CLt  = 


nr0 


(11) 


where  at  is  the  effective  diameter,  X  is  wavelength,  and  r0  is  given  by 

r0  =  3.02(C^/c2)-i  (12) 

where  C%  is  a  measure  of  turbulence  related  to  the  index  of  refraction  and  k  is  2n/X  [10]. 
At  the  range  and  frequency  of  interest,  the  effective  beam  diameter  was  found  to  be  ~5 
cm.  This  equates  to  a  turbulence  induced  jitter,  defined  as  at/R,  of  ~50  prads.  With  a 
turning  mirror  diameter  of  ~15  cm,  a  total  transmit  telescope  jitter  of  approximately  150 
prads  would  be  required  to  completely  misalign  the  system.  This  means  that  while 
turbulence  could  contribute  to  noise  in  the  baseline,  it  is  not  likely  to  cause  a  total 
misalignment  of  the  system. 
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III.  Equipment 


The  TDLAS  device  is  composed  of  the  tunable  diode  laser  source,  transmit  and 
receive  telescopes,  and  control  and  signal  collection  equipment.  A  mirror  was  included  to 
return  the  transmitted  beam  to  the  receive  telescope  so  that  the  transmit  and  receive 
telescopes  did  not  have  to  located  at  separate  sites,  easing  data  collection. 

New  Focus  Velocity  6314  and  6318  lasers  were  used  as  the  tunable  diode  sources 
depending  on  which  alkali  line  was  being  investigated,  6314  being  used  for  the  Rb  region 
and  6318  for  Cs.  Both  of  these  systems  are  capable  of  producing  a  10  mW  beam  with  a 
spectral  width  of  300  kHz.  The  6318  laser  system  has  a  full  tunable  range  of  20  nm,  from 
890  to  910  nm.  Fine  tuning  in  0.15  nm  increments  can  be  accomplished  through  the  use 
of  a  piezo  which  allows  for  fine  tuning  of  the  cavity  length.  The  6314  laser  system  has  a 
full  range  of  10  nm,  from  790  to  800  nm,  and  a  fine  tuning  range  of  approximately  0.1 
nm. 

An  Optics  for  Research  IO-5  optical  diode  was  placed  immediately  after  the  laser 
head  in  order  to  prevent  any  retroreflection  into  the  tunable  laser  cavity.  After  being 
reflected  by  a  92/8  2  pm  pellicle  beam  splitter,  a  portion  of  the  beam  was  sent  to  a  High 
Finesse  WSU-2  wavemeter.  When  used  in  combination  with  a  SIOS  SL-03  stabilized 
HeNe  laser  with  an  accuracy  of  10  MHz  for  calibration,  this  device  enabled  frequency 
measurements  accurate  to  better  than  2  MHz. 

The  remainder  of  the  beam  was  then  modulated  by  a  Thor  Labs  EO-AM-NR-C1  optical 
modulator  at  100  kHz  and  the  modulated  beam  was  coupled  to  an  RC  Optical  System 
(RCOS)  Ritchey-Chretien  12.5”  telescope  by  way  of  a  single  mode  optical  fiber. 
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Another  identical  pellicle  splitter  was  attached  to  the  telescope  between  the  output  of  the 
single  mode  fiber  and  the  telescope  itself  in  order  to  send  a  portion  of  the  beam  to  a  Thor 
Labs  PDA  100A  Si  reference  detector.  The  remainder  of  the  beam  was  expanded  to 
nearly  fill  the  secondary  mirror  of  the  telescope  through  the  use  of  an  asymmetric  lens 
with  a  focal  length  of  25  mm.  The  beam  was  then  expanded  by  the  telescope  and 
transmitted  across  a  distance  of  approximately  1  km  before  being  reflected  back  by  a  6” 
protected  gold  mirror  with  a  surface  of  flatness  of  A/20.  The  reflected  beam  was  received 
by  another  RCOS  12.5”  telescope  and  focused  onto  a  second  identical  detector,  resulting 
in  a  total  open  path  of  2,032  ±  2  m.  In  some  experiments,  a  Triad  Technologies  reference 
cell  containing  Cs  or  Rb  with  no  buffer  gas  was  placed  in  a  heater  block  between  the 
receive  telescope  and  the  receive  detector.  This  enabled  easy  spectral  location  of  the 
alkali  Di  doublet  in  the  processed  spectrum. 

Both  detectors  were  coupled  to  Stanford  Research  Systems  SR-850  lock-in 
amplifiers  capable  of  processing  a  signal  modulated  at  up  to  100  kHz  and  producing  a 
512  Hz  processed  signal.  Lock-in  time  constants  were  usually  set  at  30  or  100  ms  with  a 
filter  dropoff  of  around  6  dB/decade. 
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IV.  Experiment 


Physical  Setup 

Experimental  equipment  was  set  up  in  two  locations  separated  from  each  other  by 
1,016  ±  1  m  over  the  Upper  Machodoc  River  as  seen  in  Figure  2  below. 


Figure  2:  Diagram  of  the  test  site  showing  the  TDLAS  main  site  to  the 
north  and  the  mirror  site  to  the  south,  separated  by  1016  m.  Exact 
coordinates  for  each  site  are  shown.  Machodoc  River  runs  between 
the  sites,  beginning  to  the  west  and  emptying  into  the  Potomac  River 
to  the  east. 

The  primary  site  was  located  on  the  north  bank  of  the  river  at  the  southern  tip  of 

the  Dahlgren  Naval  Surface  Warfare  Station.  This  site  contained  an  Enterprise  16’ 

commercial  box  truck  on  blocks,  an  8’  by  8’  CONtainer  EXpress(CONEX)  military 

14 


shipping  container,  and  a  10’  by  10’  tent,  along  with  the  receiver  for  the  Scintec  BLS900 
scintillometer  and  the  standalone  Davis  Vantage  Pro2  weather  station. 

The  box  truck  primarily  functioned  as  a  shelter  and  staging  area  for  the  Tunable 
Diode  Laser  Absorption  Spectroscopy  (TDLAS)  device  and  testing  equipment.  The  RC 
Optical  System  Ritchey-Chretien  12.5”  transmit  and  receive  telescopes  were  placed  on 
QuickSet  QPT130  gimbals  and  tripods  located  just  inside  the  rear  opening  of  the  truck, 
limiting  any  obstructions  the  beam  would  encounter.  The  tunable  diode  laser  source  and 
TDLAS  control  and  signal  collection  equipment  were  placed  on  a  cart  deeper  inside  the 
truck.  The  truck  was  placed  on  blocks  in  order  to  minimize  platform  movement  allowed 
by  the  suspension  and  tires,  isolating  the  test  equipment  from  as  much  of  the  vibration 
and  movement  caused  by  wind  and  human  movement  as  possible.  By  doing  this,  the 
primary  sources  of  platform  jitter  were  limited  to  wind  buffeting  the  telescopes  and 
platform  flex. 

The  CONEX,  placed  next  to  the  truck,  was  originally  intended  for  use  as  a 

temporary  shelter  for  the  TDLAS  test  equipment,  however,  it  was  instead  used  as  a 

storage  shelter  for  extra  materials  due  to  its  small  size  and  flexibility  in  the  floor. 

The  10’  by  10’  tent  was  used  as  a  remote  shelter  from  which  the  telescopes  could  be 

aligned  and  the  TDLAS  equipment  could  be  controlled  without  disturbing  the  box  truck 

and  the  equipment  inside.  While  the  floor  of  the  truck  was  sturdy,  it  was  determined  that 

the  movement  of  a  person  from  one  end  of  the  truck  to  the  other  or  a  person  exiting  the 

truck  was  sufficient  to  misalign  the  laser  over  a  space  of  two  kilometers.  In  addition,  the 

tent  housed  the  Kanomax  handheld  particle  counter  and  its  control  laptop.  This  device 

was  placed  in  the  tent  due  to  the  need  to  routinely  recharge  the  alcohol  cartridge,  an 
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action  which  could  have  disturbed  the  truck  had  it  been  placed  with  the  rest  of  the 
equipment.  The  sample  collection  inlet  was  pointed  toward  the  entrance  to  the  tent, 
allowing  for  maximum  exposure  to  outside  air.  A  diagram  and  picture  of  the  main  test 
site  are  found  below  in  Figures  3a  and  3b. 
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Figure  3a:  Diagram  of  the  TDLAS  main  site.  The  Dahlgren  Naval 
Surface  Warfare  Center  is  located  to  the  north  and  the  river  is  to  the 
south. 
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Figure  3b:  Picture  of  the  TDLAS  main  site  showing  the  physical  setup 
of  the  test  equipment  at  the  main  test  site.  The  control  tent  is  on  the 
left,  the  CONEX  storage  container  is  in  the  middle,  and  the  box  truck 
containing  the  majority  of  the  test  equipment  is  on  the  right.  The  BLS 
900  Scintillometer  receiver  and  processing  unit  are  located  in  front  of 
the  CONEX. 

The  second  site  was  located  on  a  dock  on  the  south  side  of  the  river  on  the  north 
end  of  the  Churchill  explosives  area.  Here,  a  second  10’  by  10’  tent  was  used  to  shelter 
the  turning  mirror  as  well  as  the  transmitter  for  the  BLS900  scintillometer.  By  using 
permanent  power  located  at  the  dock,  it  was  possible  to  run  the  BLS  transmitter  for  24 
hours  a  day  for  the  duration  of  the  test.  In  addition,  a  large  2’  by  4’  alignment  target 
made  with  a  large  piece  of  rigid  foam  and  reflective  tape  was  placed  on  the  dock  in  order 
to  aid  in  location  of  the  transmitted  beam  from  the  primary  site  one  kilometer  away.  A 
diagram  and  a  picture  of  the  mirror  site  are  shown  in  Figures  4a  and  4b  below. 
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Figure  4a:  Diagram  of  the  mirror  dock  site.  The  dock  extends  into  the 
river  from  the  south 
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Figure  4b:  A  picture  of  the  mirror  site  showing  the  physical  layout  of 
the  test  equipment.  On  the  left  is  the  BLS  900  Scintillometer 
transmitter,  in  the  center  is  the  6”  turning  mirror,  and  on  the  right  is 
a  crude  retro  reflector  array  borrowed  from  Dahlgren  NSWC.  A 
large  foam  rectangle  with  reflective  strips  was  taped  to  the  shed 
behind  the  test  equipment  to  aid  in  alignment. 

Experimental  Procedure 

Before  an  experiment  was  run,  it  was  important  to  ensure  proper  alignment  of  the 
TDLAS  optical  system,  the  telescopes,  and  the  turning  mirror.  First,  the  power  being 
coupled  into  the  single  mode  fiber  going  to  the  telescope  was  maximized.  This  was 
accomplished  by  adjusting  the  location  of  the  objective  lens  and  the  placement  of  the 
fiber  relative  to  the  lens  until  the  maximum  laser  power  possible  was  coupled  to  the 
single  mode  fiber. 
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Second,  the  placement  of  the  single  mode  fiber  coupled  to  the  wavemeter  was 
adjusted  until  the  power  reaching  the  wavemeter  was  allowed  for  an  integration  time  of 
around  10  ms.  This  ensured  that  the  laser  power  reaching  the  fiber  was  to  be  strong 
enough  for  the  wavemeter  to  easily  process,  yet  not  so  strong  as  to  overwhelm  the 
wavemeter  detector. 

In  order  to  ease  the  alignment  of  the  transmitted  IR  beam,  a  5  mW  visible 
alignment  HeNe  was  mounted  to  the  top  of  the  transmit  telescope.  This  was 
accomplished  following  the  first  successful  alignment  of  the  transmit  and  receive 
telescopes.  Once  the  maximum  possible  signal  was  being  reflected  back  to  the  receive 
telescope,  the  visible  alignment  laser  was  turned  on  and  pointed  at  the  center  of  the 
turning  mirror.  When  the  coalignment  process  was  complete,  the  visible  laser  was 
reflected  back  towards  the  receive  telescope,  appearing  4  to  5  inches  below  the  telescope 
aperture.  This  offset  was  due  to  the  fact  that  the  alignment  laser  was  mounted  above  the 
transmitted  beam,  forcing  the  alignment  beam  to  return  below  the  received  beam  when 
both  telescopes  were  perfectly  aligned. 

Coarse  telescope  alignment  was  accomplished  with  an  individual  in  the  back  of 
the  box  truck.  After  coaligning  the  visible  alignment  HeNe  with  the  transmitted  TDLAS 
beam,  alignment  prior  to  an  experiment  primarily  consisted  of  locating  the  alignment 
beam  and  translating  the  telescope  such  that  it  hit  the  mirror.  At  this  point,  the  alignment 
beam  was  turned  off  and  the  transmit  telescope  was  finely  tuned  to  ensure  maximum 
atmospheric  transmission  of  the  reflected  beam.  This  was  accomplished  by  viewing  the 
mirror  through  the  receive  telescope  with  an  IR  viewfinder  and  finely  adjusting  the  angle 
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of  the  transmitting  telescope  until  the  intensity  returning  to  the  receive  telescope  reached 
a  maximum. 

Coarse  alignment  of  the  receive  telescope  was  accomplished  by  placing  an 
opaque  alignment  disc  at  the  image  focus  of  the  receive  telescope  and  changing  the  angle 
of  the  telescope  until  the  image  of  the  mirror  was  centered  on  the  disc.  The  disc  was  then 
replaced  with  a  Triad  Technologies  reference  cell  containing  the  alkali  of  interest  and 
detector  and  the  receive  telescope  was  finely  adjusted  to  maximize  the  signal  on  the 
receive  detector.  The  reference  cell  was  heated  to  just  above  the  melting  point  of  the 
alkali  of  interest  using  cartridge  heaters. 

A  similar  method  was  also  used  to  ensure  proper  alignment  of  the  mirror.  If 
mirror  alignment  was  required,  the  angle  of  the  mirror  was  changed  until  the  image  of  the 
transmitting  telescope  aperture  as  viewed  through  the  receive  telescope  was  centered  on 
the  mirror. 

Fine  telescope  alignment  was  then  accomplished  from  the  remote  control  tent. 
Because  of  the  slight  change  in  the  vertical  pointing  angle  caused  by  the  weight  of  an 
individual  in  the  truck,  the  transmit  and  receive  telescopes  were  adjusted  to  point  down 
slightly  after  the  truck  was  empty,  enabling  the  maximum  signal  to  be  collected  by  the 
receive  telescope. 

Atmospheric  Absorption  Spectra 

From  17  October  2012  to  23  October  2012,  full  spectra  were  collected  near  the 

Rb  and  Cs  Di  emission  lines.  In  the  case  of  Cs,  a  New  Focus  Velocity  6318  was  used  as 

the  tunable  diode  laser  source.  Two  spectral  ranges  were  investigated,  one  comprised  of 

the  entire  tunable  range  of  the  diode  and  the  other  focusing  on  the  range  5  nm  on  either 
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side  of  the  Cs  Di  line.  The  wider  range  was  typically  conducted  using  a  scan  speed  of 
approximately  0.005  nm  per  second  with  one  or  two  averages  per  piezo  range.  These 
parameters  were  selected  in  order  to  keep  collection  times  to  below  two  hours.  This  time 
limit  was  imposed  due  to  the  difficulty  in  maintaining  telescope  alignment  over  2  km  for 
an  extended  period  of  time.  Typically,  telescopes  stayed  aligned  for  between  100  to  150 
minutes.  For  the  narrower  spectral  range,  it  was  possible  to  fit  the  same  time  constraint 
while  using  a  scan  speed  as  low  as  0.003  nm  per  second  as  well  as  a  greater  number  of 
averages. 

In  the  case  of  Rb,  a  New  Focus  Velocity  6314  was  used  as  the  laser  source.  This 
system  has  a  full  scanning  range  of  from  790  to  802nm  and  a  piezo  scan  range  of  about 
0.15  nm.  Because  of  the  reduced  full  tunable  range,  the  narrow  range  was  selected  to  be 
only  2.5  nm  on  either  side  of  the  Rb  Di  line.  Scan  rates  were  the  same  as  those  used  in 
the  Cs  case. 

Long  Term  Single  Frequency  Baseline  Monitoring 

In  order  to  quantify  how  the  standard  deviation  of  the  baseline  changed  as  a 
function  of  turbulence  or  scattering,  a  single  frequency  located  away  from  atmospheric 
absorption  features  was  monitored  continuously  for  a  period  of  one  or  two  hours.  Due  to 
the  lack  of  strong  FLO  absorption  features  in  its  scan  range,  the  New  Focus  6314  diode 
laser  was  used  in  this  experiment.  The  frequency  selected  for  monitoring  was  12599.34 
cm"1.  Due  to  the  long  term  nature  of  the  experiment,  data  was  saved  to  the  hard  drive 
every  ten  minutes  to  prevent  the  total  loss  of  data  during  an  experiment  due  to  a  system 
malfunctions  such  as  software  crashes  or  hardware  misalignments. 
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Long  Term  Single  Line  Monitoring 

In  order  to  extract  temperature  and  concentration  as  a  function  of  time,  two 
moderately  strong  rotational  FLO  absorption  lines  were  monitored  over  the  course  of  an 
hour.  These  lines  were  chosen  based  on  their  strength  and  resolvability  from  other 
nearby  lines.  In  order  to  avoid  saturating  the  lock-ins,  which  can  only  reliably  detect 
absorption  to  around  4,  the  selected  lines  were  no  greater  than  a  2  in  absorbance.  Based 
on  these  criteria,  the  lines  located  at  11192.69  and  11176.80  cm  1  were  selected  for 
monitoring.  The  New  Focus  Velocity  6318  laser  head  was  chosen  for  this  experiment 
due  to  the  presence  of  several  strong  FLO  lines  in  its  operating  range.  In  order  to  obtain 
data  for  both  lines  as  near  to  each  other  in  time  as  possible,  each  line  was  scanned  for  one 
minute,  after  which  the  laser  would  tune  to  the  second  line,  scanning  it  for  one  minute 
before  returning  to  the  first  line.  This  process  was  repeated  up  to  30  times,  gathering 
pairs  of  lines  that  could  be  fit  LBLRTM  spectra  generated  using  the  HITRAN  2008 
database  in  order  to  extract  temperature  and  concentration  estimates  that  could  be 
compared  with  collected  weather  data 
Meteorological  Data 

Turbulence  data  was  collected  by  a  BLS900  Boundary  Layer  Scintillometer.  The 
transmitter,  placed  under  the  tent  at  the  mirror  site,  was  mounted  on  a  tripod  1.2  m  above 
the  ground.  In  order  to  reliably  obtain  the  maximum  range  of  crosswind  speeds,  0. 1  to  30 
m/s,  the  maximum  possible  pulse  repetition  rate  of  125Hz  was  used.  Permanent  power  at 
the  mirror  site  allowed  for  collections  to  be  taken  for  the  entire  duration  of  the  test.  The 
receiver,  located  at  the  main  TDLAS  site,  was  also  placed  on  a  tripod  1.2  m  off  the 

ground.  The  signal  processing  unit  was  placed  on  the  ground  next  to  the  receiver  tripod. 
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At  1016  m,  the  scintillometer  was  able  to  measure  Cn  turbulence  values  between 
4  x  10  15  and  4  x  1 0  1 1 .  The  integration  time  of  the  processing  unit  was  set  to  the 
minimum  allowed  value  of  1  minute  in  order  to  maximize  the  number  of  data  points 
collected. 

Particulate  information  was  collected  by  a  Kanomax  handheld  particle  counter 
placed  in  the  remote  control  tent.  This  device  is  capable  of  measuring  particles 
concentrations  of  between  0  and  100,000  particles/cm  for  particles  ranging  in  size  from 
0.015  to  1  pm.  This  device  was  set  to  collect  data  once  every  five  seconds  in  four  hour 
spans.  Typical  particle  counts  were  between  2,000  and  10,000  particles/cm3  with  one 
instance  of  counts  above  50,000  particles/cm  corresponding  to  user  actions  around  the 
counter. 

Weather  data  was  collected  by  a  Davis  Vantage  Pro2  Weather  Station  set  up  at  the 
main  site.  This  station  consists  of  a  console  unit  placed  inside  the  truck  and  a  remote  unit 
outside.  This  device  collected  continuous  data  for  relative  humidity,  pressure,  and 
temperature  at  15  minute  intervals  with  95%  confidence  bounds  of  6%,  2  hPa,  and  1  C 
respectively. 

Testing  Schedule 

15  October,  Set  Up,  no  collection  (Started  BLS  Data  Collection) 

16  October,  1600-2100  -  Cs  Long  Term  Absorption  Monitoring 

17  October,  1500-2300  -  Cs  Spectra 

18  October,  1100-1230  -  Cs  Spectra 

19  October,  1300-2359  -  Rb  Spectra 
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20  October,  0000-0100  -  Rb  Spectra  (continued  from  19  Oct) 

20/21  October,  Data  Analysis  Days,  No  Collection 
22  October,  1500-1830  -  Rb  Spectra 

22  October,  1900-2100  -  Cs  Spectra 

23  October,  1330-1700  -  Rb  Spectra 

23  October,  1730-1930  -  Rb  Long  Term  Baseline 

24  October,  1 130-1930  -  Rb  Long  Term  Baseline 

25  October,  Pack,  No  Collection  (Ended  BLS  Data  Collection) 

Table  1:  Schedule  showing  what  data  was  collected  on  each  day  of  the  test.  Data 
was  collected  between  16  October  and  25  October,  with  a  day  before  and  after  this 
for  set  up  and  packing  respectively.  Cs  and  Rb  Full  represent  days  where  full 
spectra  were  collected  for  each  species.  Cs  Long  represents  a  long  term  absorption 
line  measurement  and  Rb  Long  represents  a  long  term  baseline  monitoring 
collection. 


Sun _ Mon _ Tuesday _ Wed _ Thur _ Fri _ Saturday 


14  Oct 

15  Oct 

16  Oct 

17  Oct 

18  Oct 

19  Oct 

20  Oct 

Travel 

Day 

Set  Up 
Day 

1600-2100 
Cs  Fong 

1500-2300 
Cs  Full 

1100-1230 
Cs  Full 

1300-2359 
Rb  Full 

0000-0130 
Rb  Full 

21  Oct 

22  Oct 

23  Oct 

24  Oct 

25  Oct 

26  Oct 

27  Oct 

Data 

Day 

1500-1830 
Rb  Full 
1900-2100 
Cs  Full 

1330-1700 
Rb  Full 
1730-1930 
Rb  Fong 

1130-1930 
Rb  Fong 

Pack 

Day 

Travel 

Day 
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Y.  Results 


Meteorological  Data 

A  BLS  900  Scintillometer  was  used  to  collect  turbulence  measurements  for  the  9 
day  duration  of  the  experiment.  At  a  range  of  1,000  meters,  the  detectable  range  of  Cn~ 
values  was  5  x  10  15  to  10  12.  A  plot  of  the  turbulence  data  for  all  9  days  as  well  as  a 
sample  from  a  single  day  are  found  in  Figures  5  a  and  5b. 


Figure  5a:  Cn2  measurements  collected  over  the  duration  of  the 
experiment.  The  turbulence  values  cycle  from  ~5  x  10'13  at  night  to  ~5 
x  10  ls  in  the  afternoon  just  prior  to  sunset. 
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Figure  5b:  Turbulence  measurements  for  21  October,  highlighting  the 
drastic  rise  in  turbulence  immediately  following  sunset. 

15  12 

Observed  turbulence  measurements  ranged  from  -5x10"  to  10"  ,  within  the 
detectable  range  of  the  BLS  900  Scintillometer  at  1,000  meters.  The  peak  typically 
occurred  in  mid-morning,  with  Cn~  values  dropping  through  the  afternoon  to  a  low  of 
approximately  10  14  just  before  sunset.  The  turbulence  values  consistently  spiked  by  an 
order  of  magnitude  immediately  after  sunset  before  slowly  rising  through  the  night  and 
early  morning  to  peak  again  in  mid-morning.  While  turbulence  on  the  order  of  that 
measured  is  not  sufficient  to  misalign  or  completely  degrade  the  propagating  signal,  it  is 
possible  that  data  taken  during  more  turbulent  periods  could  have  a  noisier  baseline. 
Particulate  counts  were  collected  for  several  four  hour  time  spans  over  the  course  of  the 
experiment.  The  four  hour  limit  is  dictated  by  the  length  of  time  a  single  alcohol 
cartridge  would  function  before  having  to  be  recharged.  A  sample  of  two  of  these  scans 
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taken  in  quick  succession  on  the  afternoon  of  22  October  is  shown  in  Figure  6  below. 


Time,  22  October, 

Figure  6:  Particle  counts  for  a  seven  hour  block  of  time  on  22  October 
2012.  This  sample  includes  a  count  fall  off  around  17:15 
corresponding  to  the  alcohol  cartridge  in  the  particle  counter  failing 
as  well  as  a  spike  of  over  50,000  particles/cm3  corresponding  to  the 
spray  of  insect  repellant  near  the  collector  inlet.  Typical  particle 
counts  were  on  the  order  of  1,000  to  2,000  particles/cm3. 

Typical  particle  values  for  the  duration  of  the  experiment  were  on  the  order  of 

-3 

10,000  to  20,000  particles/cm  .  The  maximum  possible  collection  time  was  ~4  hours,  at 
which  point  the  alcohol  cartridge  must  be  removed  from  the  counter  and  recharged  for 
half  an  hour.  The  maximum  value  recorded  was  over  50,000  particles/cm  ,  which 
corresponds  to  a  rapid  influx  of  insect  spray  near  the  inlet  of  the  particle  counter.  Several 
small  peaks  such  as  the  one  shortly  after  1800  were  observed  and  attributed  to  vehicle 
traffic  around  the  test  site. 
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Atmospheric  Absorption  Spectra 

Broad  Cesium  Region  Spectra 

Data  were  collected  for  two  spectral  ranges,  one  consisting  of  the  entire  tunable 
range  of  the  laser  diode  and  another  focused  on  the  range  around  the  DPAL  Di  emission 
line.  In  the  Cs  emission  region,  the  maximum  tunable  range  was  from  890.5  to  909.9  nm, 
dictated  by  the  tunable  range  of  the  laser  diode  used.  An  example  of  the  spectrum 
collected  over  this  range  is  found  in  Figure  7  below. 


11,000  11,020  11,040  11,060  11,080  11,100  11,120  11,140  11,160  11,180  11,200  11,220 

Wavenumber,  u,  (cm-1) 


Figure  7:  Atmospheric  absorption  spectrum  of  the  entire  20  nm  range 
of  the  New  Focus  Velocity  6318  Tunable  Diode  Laser.  The  majority 
of  the  spectrum  is  made  up  of  the  rotational  structure  of  the  (0  0 
0) — >(0  0  3)  vibrational  transition  with  some  small  contributions  by  the 
(000) — >(  102)  transition. 

This  spectrum  was  collected  at  night,  from  8:07  PM  to  11:13  PM,  on  18  October 
2012  and  is  dominated  by  the  R-branch  of  the  (0  0  0)  — >  (0  0  3)  vibrational  band  of  H2O 
with  rotational  quantum  numbers  ranging  from  0  to  10,  the  majority  of  which  exhibit  an 
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Absorption  Coefficient,  a(v),  (km  1 


absorbance  of  2  or  greater  at  this  distance.  The  P-branch  and  Q-branch  have  significant 
contributions  at  the  lower  energy  end  of  the  scan.  Several  weaker  contributions  by  the  (0 
0  0)  — >  ( 1  0  2)  vibrational  band  of  H2O  are  scattered  throughout  the  entire  scan,  all  with 
an  absorbance  of  less  than  0.5.  In  total,  over  250  absorption  lines  were  identified  in  this 
spectrum.  The  locations  of  these  lines,  along  with  their  relevant  HITRAN  data  are  listed 
in  Table  4  of  Appendix  A. 

The  spectrum  was  fit  to  the  known  spectral  lines  in  the  HITRAN  database  using 
LBLRTM  by  varying  the  temperature,  pressure,  and  concentration  estimates  until  the 
residuals  of  a  Voigt  fit  were  minimized.  The  result  was  a  temperature  estimate  of  282.6  ± 
0.1  K,  a  pressure  estimate  of  0.9590  ±  0.0004  atm,  and  a  concentration  estimate  of  3.2350 
±  0.0009  x  10  cm"  ,  with  95%  confidence  intervals  on  all  three  uncertainties.  The 
standard  deviation  of  the  residuals  of  this  final  model  fit  is  a  measure  of  the  noise  in  the 
baseline  and  determines  the  minimum  detectable  absorbance. 

With  a  baseline  noise  of  0.14  and  a  maximum  detected  absorbance  of  1 1.7,  the 
signal  to  noise  ratio  (SNR)  of  this  data  set  was  calculated  to  be  83.6.  Due  to  saturation  of 
the  lock-ins,  this  SNR  is  not  entirely  accurate.  If  instead  the  maximum  detectable 
absorbance,  4,  is  used,  the  SNR  becomes  28.6.  With  a  maximum  observed  turbulence 
value  of  3  x  10"14,  the  turbulence  induced  jitter  in  this  case  was  calculated  to  be  16  prads. 
This  is  relatively  low  compared  to  the  jitter  induced  by  wind,  which  averaged  2.98  m/s 
over  the  course  of  this  collection.  While  neither  wind  at  this  speed  nor  turbulence  is 
sufficient  to  permanently  misalign  the  system,  the  jitter  induced  by  buffeting  of  the  truck 
itself,  the  telescopes  in  the  truck,  and  the  folding  mirror  by  the  wind  can  induce 

significant  noise  levels  into  the  data  if  it  is  not  mitigated  by  signal  averaging. 
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Narrowed  Cesium  Region 


In  one  experiment,  on  22  October  2012  at  7:05,  a  cell  containing  cesium  at  34°  C 
was  included  in  the  open  path  of  the  laser.  This  was  done  in  order  to  easily  identify  the 
precise  location  of  the  cesium  Di  emission  line  at  894.6  nm.  The  full  range  of  this  scan 
was  from  891.7  to  897.5  nm,  approximately  3  nm  on  either  side  of  the  Di  line.  The  3.5 
cm  1  portion  of  this  scan  immediately  surround  the  Di  line  is  shown  in  Figure  8. 


Wavenumber,  />  —  vref,  (cm  *) 


Figure  8:  0.3  nm  sample  of  the  atmospheric  absorption  spectrum 
collected  over  the  6  nm  range  immediately  surrounding  the  Cs  Di 
doublet.  vref  is  the  energy  of  the  lowest  frequency  line  of  the  Cs 
doublet,  11178.112  cm'1.  The  dashed  line  represents  the  LBLRTM  fit 
to  the  collected  data,  displayed  as  black  points.  The  Cs  doublet 
appears  as  a  dotted  line  and  was  excluded  from  the  atmospheric 
model  fit  as  Cs  parameters  are  not  included  in  atmospheric  databases. 
The  gaps  in  the  spectrum  are  a  product  of  the  tunable  laser  failing  to 
completely  overlap  consecutive  scans. 

The  gray  data  points  correspond  to  the  Cs  Dj  doublet.  This  data  was  excluded 

from  the  nonlinear  fit  because  the  HITRAN  database  does  not  include  the  Cs  spectral 
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lines.  With  a  standard  deviation  in  the  residual  of  0.1 1  and  a  maximum  detected 
absorbance  of  8.4,  the  signal  to  noise  ratio  of  this  entire  scan  was  calculated  to  be  76.4. 
Two  FTO  absorption  features  are  clearly  visible  on  either  side  of  the  Di  doublet,  both 
corresponding  to  the  rotational  structure  of  the  (0  0  0)  — >  (0  0  3)  vibrational  transition. 
One  line,  corresponding  to  the  [8  4  5]  — >  [9  4  6]  rotational  transition,  is  located  at  894.7 
nm.  The  second  feature  is  made  up  of  two  lines,  one  of  which  is  due  to  the  [9  5  5]  — >  [10 
5  6]  rotational  transition  and  the  other  of  which  is  due  to  [6  2  4]  — >  [7  2  5]  rotational 
transition  and  is  located  at  894.45  nm. 

Fit  Correlations  with  Meteorological  Data 

Full  absorption  spectra  were  recorded  throughout  the  9  days  of  data  collection, 
resulting  in  6data  sets  from  which  fit  estimates  for  temperature,  pressure  and 
concentration  could  be  extracted.  In  addition,  the  Davis  Vantage  Pro2  weather  station 
recorded  temperature  and  pressure  for  the  duration  of  the  test.  The  temperature,  pressure, 
dew  point,  and  relative  humidity  provided  by  the  Davis  weather  station  were  used  to 
obtain  concentration  measurements  for  the  duration  of  the  test.  Figures  9a  and  9b  below 
show  how  the  fit  estimates  compare  to  the  measured  temperature  and  pressure  values 
respectively. 
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H-O  Concentration,  (1017  cm  3)^  Temperature,  (C) 


Figure  9a:  Temperature  measurements  from  the  Davis  Vantage  Pro2 
weather  station  (dark  grey  line),  national  weather  station  in  Quantico, 
VA  (dashed  line),  and  spectral  fit  estimates  (data  markers)  for  the  9 
day  duration  of  data  collection. 


Figure  9b:  Temperature  measurements  from  the  Davis  Vantage  Pro2 
weather  station  (dark  grey  line),  national  weather  station  in  Quantico, 
VA  (dashed  line),  and  spectral  fit  estimates  (data  markers)  for  the  9 
day  duration  of  data  collection. 
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The  data  sets  used  in  these  fit  estimates  come  from  two  collections  on  the  evening 
of  17  October,  two  on  the  evening  of  19  October,  one  on  the  evening  of  22  October,  and 
one  from  23  October.  Because  of  the  broad  spectral  region  scanned  on  17  October,  each 
of  these  scans  was  broken  up  into  3  parts  and  each  region  was  fitted,  providing  three 
more  data  points  per  scan.  This  process  was  also  applied  at  the  test  on  the  23rd. 

The  temperature  estimates  from  the  experimental  fits  vary  more  widely  and  faster 
than  the  Davis  weather  station  measurements  with  estimates  in  the  evening  of  18  October 
dropping  at  a  higher  rate,  3.18  degrees  per  hour,  and  ending  at  a  lower  temperature,  6.6° 
C,  than  the  weather  measurements  at  0.73  degrees/hour  and  11.9°  C  respectively. 

Overall,  temperature  estimates  ranged  from  6.60°  to  30.03°  C  while  the  Davis 
weather  station  reported  values  between  6.45°  and  26.75°  C,  with  an  agreement  of  better 
than  6%. 

In  two  cases,  it  was  possible  to  divide  the  collected  spectrum  into  3  pieces  which  were 
then  fit  to  LBLRTM  independently.  One  of  these  scans,  taken  on  17  October,  took  place 
over  a  steep  drop  in  temperature,  allowing  for  a  comparison  estimates  with  the  Davis 
weather  station  on  a  faster  time  scale.  The  estimates  and  errors  for  these  three  pieces  as 
well  as  the  entire  scan  are  shown  in  Figure  10  below. 
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Figure  10:  Fit  estimates  and  weather  station  data  for  temperature  on 
17  October  2012.  Vertical  error  bars  correspond  to  95%  statistical 
errors  from  the  fits  and  horizontal  error  bars  represent  the  period  of 
time  covered  by  each  piece  of  the  collection.  The  solid  grey  line 
represents  the  Davis  weather  station  data  and  the  dashed  line 
represents  the  temperature  data  from  the  national  weather  station  at 
Quantico. 

Concentration  estimates  for  all  of  the  analyzed  scans  ranged  from  2.55  x  10  and 

17  3 

3.5  x  10  cm"  with  individual  values  being  between  1  and  25%  lower  than  calculated 
values  which  ranged  from  1.79  x  10  to  5.32  x  10  cm'  he  estimates  which  were 
furthest  from  the  calculated  values  were  largest  at  times  when  the  concentration  was 
changing  most  rapidly.  This  discrepancy  is  likely  due  to  the  Davis  weather  station  being 
placed  several  meters  away  from  the  edge  of  the  river.  Because  the  TDLAS  laser  is 
interrogating  air  over  the  river,  it  is  possible  that  the  both  the  absolute  concentration  and 
the  rate  at  which  that  concentration  changes  as  measured  by  the  TDLAS  device  could 
differ  from  the  value  calculated  by  the  meteorological  data.  In  times  of  more  stable 
concentrations,  the  discrepancy  was  reduced  to  between  1  and  10%.  A  summary  of  the 
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temperature  and  concentration  estimates  for  each  of  the  scans  plotted  is  found  in  Table  2 
below.  In  addition,  plots  of  all  data  sets  used  to  obtain  these  estimates  are  found  in 
Appendix  B. 


Table  2:  Scan  times,  fit  estimates,  95%  confidence  bounds  (+),  Davis  weather 
station,  and  percent  error  for  each  scan  analyzed. 

Time  Concentration  (1017  cm'3)  Temperature  (K) 


Start  End  Fit  ±  Davis  Error  Fit  +  Davis  Error 


Comparison  with  Previous  Work 

Recent  research  has  been  performed  using  the  same  TDLAS  device  used  in  this 
experiment  to  investigate  Cs  and  Rb  emission  regions  [11-15].  These  experiments  were 
performed  at  150  m  indoors  and  at  1  km  outdoors  on  the  Air  Force  Institute  of 
Technology  (AFIT)  Laser  Range.  This  previous  work  was  able  to  use  extensive  signal 
averaging  to  drive  the  signal  to  noise  ratio  higher  than  1,000.  Such  a  high  ratio  is 
important  in  the  accurate  estimation  of  temperature,  pressure,  and  concentration.  Figure 
1 1  below  shows  a  spectrum  collected  at  2,000  m  along  with  these  previous  spectra, 
collected  at  150  and  1,000  m. 
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Figure  11:  Atmospheric  absorption  spectra  collected  over  the  Cs  Di 
spectral  region  at  150  m  (black),  1,000  m  (dark  gray),  and  2,000  m 
(light  gray).  The  D|  doublet  is  included  in  the  spectrum  collected  at 
150  m.  The  increase  in  baseline  noise  is  clear  as  the  experimental  path 
length  is  increased,  with  a  dramatic  rise  going  to  2,000  m  due  to  a  lack 
of  signal  averaging  as  well  as  additional  platform  jitter. 

The  noise  in  the  spectrum  collected  at  the  2,000  meter  path  length  is  significantly 

higher  than  that  at  either  previous  distance.  There  are  several  reasons  for  this 

degradation.  First,  the  wind  speeds  for  this  collection  were  around  3  m/s.  This  can 

induce  jitter  in  the  mirror  as  well  as  in  the  platform  the  telescopes  were  on.  At  around  50 

prads,  turbulence  jitter  could  have  also  contributed  to  the  degradation  of  the  baseline, 

introducing  noise  without  misaligning  the  system. 

Further,  by  observing  the  intensity  of  the  line  at  11182.5  cm  1  at  each  range,  1.73, 

0.44,  and  0.12  at  2  km,  1  km,  and  150  m  respectively,  it  is  seen  that  the  strength  of  the 

line  at  2  km  is  3.9  times  that  of  the  1  km  distance.  In  addition,  the  intensity  at  2  km  is 
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13.3  times  that  at  the  150  m  case.  Because  of  Beer’s  Law,  it  is  known  that  absorbance  is 
linearly  related  with  both  path  length  and  concentration.  Because  absorption  varies 
linearly  with  concentration,  it  is  likely  that  the  concentration  of  water  molecules  at 
Dahlgren,  VA  in  October  is  roughly  twice  that  at  AFIT  in  winter,  the  location  and  date  of 
the  1  km  experiment.  The  intensity  ratio  between  the  150  m  and  2  km  lines,  however,  is 
only  slightly  higher  than  the  ratio  of  the  distances,  meaning  the  FLO  concentrations  were 
likely  similar  during  those  tests. 

These  observations  are  consistent  with  concentration  measurements  taken  at  all 

IV 

three  test  sites.  For  the  150  m  path  length,  a  concentration  of  3.096  ±  0.004  x  10 

3  17  3 

molecules/cm  was  recorded.  A  concentration  of  1.94  ±  0.004  x  10  molecules/cm  was 

17 

reported  for  the  1  km  path  length.  A  concentration  estimate  of  3.23  ±  0.001  x  10 
molecules/cm  was  extracted  from  the  data  at  the  time  of  the  2  km  test.  The  small 
increase  in  concentration  at  the  2  km  test  with  respect  to  the  150  m  test  makes  up  the 
small  discrepancy  seen  in  the  intensity  ratio.  The  distance  ratio  multiplied  by  the 
concentration  ratio  between  the  two  tests  returns  a  total  expected  ratio  of  13.9,  only  3.4% 
below  the  observed  value.  In  the  case  of  the  1  km  test,  the  total  calculated  ratio  was  3.33, 
15%  below  the  observed  intensity  ratio.  The  remaining  discrepancies  in  intensity  are 
likely  due  to  other  meteorological  parameters,  specifically  temperature.  The  temperature 
reported  by  the  Davis  weather  station  reported  a  value  of  16.9  C,  slightly  higher  than  the 
15  C  reported  for  the  150  m  case  and  significantly  higher  than  the  6  C  reported  for  the  1 
km  case.  A  summary  of  the  concentration,  path  length,  and  absorption  ratios  can  be 
found  in  Table  3  below. 
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Table  3:  Observed  concentration  and  absorption  values  for  150,  1,000,  and  2,000 
meter  collections.  All  ratios  are  with  respect  to  the  relevant  values  at  2,000  meters. 
The  ratio  product  is  determined  by  multiplying  the  path  ratio  by  the  concentration 
ratio. 


Length 

(m) 

[H2o] 

(x  1017  cm'3) 

Observed 

Absorption 

Path 

Ratio 

[H2o] 

Ratio 

Ratio 

Product 

Intensity 

Ratio 

Error 

150 

3.096 

0.12 

0.075 

0.932 

0.0699 

0.0694 

0.73 

1,000 

1.940 

0.44 

0.500 

0.601 

0.3003 

0.2543 

18.07 

3.230 

1.73 

1 

1 

1 

1 

0 

Rotational  Dependence  of  Broadening  Rates 

HITRAN  lists  values  for  air  and  self-broadened  half  widths  as  well  as  a  method 
by  which  to  correct  these  values  to  report  an  actual  observed  half  width  [16].  There  is  a 
slight  change  in  the  broadening  coefficient  that  corresponds  with  a  change  in  the 
rotational  transition  being  observed.  By  finding  two  widely  separated  rotational  lines  in 
the  same  vibrational  transition,  it  was  possible  to  investigate  whether  or  not  it  was 
possible  to  observe  this  change  in  the  data  collected  at  Dahlgren.  The  two  lines  selected 
were  corrected  to  be  centered  at  zero  wavenumbers  and  normalized  such  that  the  peak 
absorbance  was  1,  as  seen  in  Figure  12  below. 
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1.2 


Figure  12:  Absorption  lines  corresponding  to  the  [1  1  0]  — >  [2  1  1] 

(light  gray)  and  [6  2  5]  — >  [7  2  6]  (dark  gray)  rotational  transitions 
within  the  vibrational  (0  0  0)  — ►  (0  0  3)  band,  overlaid  and  normalized 
to  a  peak  intensity  of  1.  An  LBLRTM  fit  for  each  line  is 
superimposed  on  the  data,  with  the  dotted  line  corresponding  to  the  J 
=  2  line  and  the  dashed  line  corresponding  to  the  J  =7  line.  Noise  in 
the  top  of  each  data  set  is  due  to  absorbance  beyond  the  detectable 
limit  of  the  lock-in  amplifiers. 

Based  on  this  plot,  it  was  determined  that  the  half  width  of  the  J=2  rotational  line 
was  0.1 1  cm’1.  HITRAN  reports  this  width  as  0.082  cm’1,  leaving  a  discrepancy  between 
the  reported  value  and  the  experimental  value  24%.  The  J=7  rotational  line  had  a  width 
of  0.08  cm  1  28%  higher  than  the  value  reported  in  HITRAN,  0.059  cm’1.  It  is  likely  that 
the  large  discrepancy  between  the  observed  and  theoretical  values  stems  from  the  noise  in 
the  data.  Because  the  determination  of  pressure,  and  thus  linewidth,  is  so  dependent  on 
the  wings  of  each  absorption  feature,  there  is  a  relatively  large  amount  of  uncertainty 


induced  by  noise  in  the  data. 
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HITRAN  lists  the  uncertainties  in  the  reported  values  for  the  air-broadened  half¬ 
width  as  being  between  5  and  10%.  The  values  reported  for  the  self-broadened  half 
width  are  listed  as  being  either  an  average  of  multiple  experimental  results  or  as  a  best 
estimate,  with  an  unknown  amount  of  error.  In  addition,  the  amount  of  noise  present  in 
the  experimentally  collected  2  km  data  increases  the  error  inherent  in  the  nonlinear  fit  as 
well  as  complicating  the  determination  of  the  half  width  from  the  data  itself,  with  the 
range  of  frequencies  corresponding  to  a  normalized  intensity  of  0.5  being  approximately 
0.015  cm1.  Despite  the  error  in  the  broadening  rates  themselves,  the  ratio  of  the 
broadening  rates  at  the  J=2  line  to  the  line  at  J=7  is  1.313  as  determined  by  the  fit  values, 
only  5%  below  the  ratio  of  the  values  reported  in  HITRAN. 
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VI.  Conclusions  and  Recommendations 


Conclusions  of  Research 

A  ruggedized  TDLAS  device  has  been  deployed  to  Dahlgren,  VA  to  investigate 
atmospheric  absorption  at  DPAL  wavelengths  in  a  maritime  environment  at  a  path  length 
of  2  km.  Over  the  course  of  9  days,  experiments  were  conducted  to  investigate 
absorption  around  the  DPAL  emission  lines  of  Rb  and  Cs,  the  way  in  which  changing 
water  vapor  concentrations  can  affect  the  strength  of  absorption  features,  and  the  effects 
of  various  parameters  on  the  noise  in  the  baseline  as  well  as  the  value  of  the  baseline 
itself.  This  series  of  experiments  showed  the  feasibility  of  later  investigations  at  still 
longer  path  lengths  as  well  as  highlight  some  of  the  challenges  that  would  accompany 
such  an  effort.  In  addition,  performing  this  investigation  at  a  longer  path  length  showed 
how  absorption  features  grow  and  broaden  as  the  propagation  distance  grows.  This 
information  is  valuable  in  determining  how  molecular  absorption  will  affect  a 
propagating  laser  beam. 

Spectra  were  collected  within  the  20  nm  range  surrounding  the  Cs  Di  laser 
emission  line  as  well  as  in  a  10  nm  range  around  the  Rb  D[  line.  With  a  maximum 
detectable  absorbance  of  around  4  and  a  baseline  noise  level  of  around  0.1,  the  SNR  for 
the  spectral  region  surrounding  Cs  was  limited  to  40.  With  a  lower  peak  detected 
absorbance  due  to  the  relative  strength  of  the  H2O  absorption  lines,  the  Rb  spectral  region 
typically  had  a  SNR  of  around  10.  These  values  are  significantly  lower  than  those  from 
previous  efforts  using  the  same  device,  as  high  as  1,000,  due  to  the  extended  path  length, 
increased  platform  jitter,  and  the  lack  of  signal  averaging.  Temperature  and 
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concentration  estimates  were  made  by  fitting  the  collected  spectra  to  the  HITRAN 
database  using  LBLRTM.  These  estimates  were  used  in  conjunction  with  collected 
weather  data  to  investigate  potential  differences  in  the  ways  temperature  and 
concentration  change  over  time  over  land  and  water.  Temperature  estimates,  which 
ranged  from  6.6°  C  at  night  to  30.02°  C  in  the  mid-aftemoon,  were  within  6%  of  the 
corresponding  values  reported  by  a  Davis  Vantage  Pro2  weather  station,  which  ranged 
from  6.5°  to  26.8°  C.  In  addition,  some  differences  in  the  rate  at  which  the  temperature 
changed  over  water  as  opposed  to  on  land  were  noted,  with  the  air  over  the  water 
appearing  to  cool  around  4  times  faster  than  the  air  over  the  land.  Concentration 
estimates  were  within  25%  of  the  values  calculated  from  the  collected  meteorological 
data,  with  the  largest  difference  being  found  at  points  in  time  when  the  concentration  was 
changing  most  drastically.  Concentration  estimates  ranged  from  2.55  x  10  and  3.5  x 
10  cm"  with  calculated  values  for  the  entire  9  day  period  of  data  collection  ranging 
from  1.79  x  1017  to  5.32  x  1017cm'3. 

Spectral  and  meteorological  data  from  this  2  km  investigation  were  compared  to 
data  collected  at  path  lengths  of  150  m  and  1  km.  By  overlaying  the  collected  spectra,  it 
was  possible  to  directly  observe  how  the  absorption  features  change  as  a  function  of  path 
length  and  hBO  concentration.  Because  the  water  vapor  concentrations  were  similar,  on 
the  order  of  3  x  10  molecules/cm  ,  between  the  150  m  and  the  2  km  tests,  the  relative 
strength  of  these  absorption  features  was  dictated  primarily  by  the  change  in  path  length. 
Because  the  concentration  at  the  1,000  m  test  was  roughly  half  that  at  the  other  tests,  the 
relative  strength  of  the  absorption  features  was  dictated  in  part  by  the  change  in  path 

length  and  in  part  due  to  the  change  in  species  concentration. 
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By  analyzing  the  linewidths  of  two  rotational  lines,  the  [110]  — ►  [2  11] 
transition  at  1 1080  cm  1  and  the  [6  2  5]  — >  [7  2  6]  transition  at  11153  cm"1,  from  the  (0  0 
0)  to  (0  0  3)  vibrational  band,  it  was  possible  to  investigate  the  rotational  dependence  of 
broadening  rates.  These  broadened  linewidths,  as  well  as  the  ratio  by  which  the 
broadening  rate  changes  as  a  function  of  rotational  quantum  number,  were  compared  to 
the  values  reported  in  the  HITRAN  database  for  the  given  transitions.  The  broadening 
rates  from  the  two  investigated  lines  were  found  to  be  0.109  cm  1  for  the  lower  line  and 
0.083  cm  1  for  the  upper  line,  both  of  which  are  between  25  and  30%  higher  than  those 
reported  by  the  HITRAN  database,  0.0823  and  0.0594  cm"1  respectively.  Despite  this 
discrepancy,  the  ratio  by  which  the  broadening  rate  decreases  as  a  function  of  rotational 
quantum  number  in  the  collected  data  is  within  5%  of  ratio  of  the  values  found  in  the 
database. 

Recommendations  for  Future  Research 

Continuous  upgrades  to  the  TDLAS  device  used  in  this  experiment  have  allowed 
for  a  gradual  increase  in  its  utility.  By  upgrading  the  tunable  diode  laser  system,  the 
mirror,  and  the  signal  collection  and  processing  equipment,  the  capabilities  of  the  AFIT 
TDLAS  apparatus  will  continue  to  increase  dramatically. 

In  terms  of  the  tunable  diode  laser  system  used  to  prepare  the  beam  prior  to 
expansion,  replacing  the  free  space  system  currently  used  with  a  fiber  coupled  scheme 
would  have  several  benefits.  First,  it  would  drastically  cut  down  on  the  amount  of  time 
required  to  change  diode  laser  heads,  increasing  the  efficiency  of  future  experiments.  In 
its  current  configuration,  the  changing  of  a  laser  head  alone  takes  15-30  minutes,  with  a 
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full  alignment  of  all  optics  taking  several  hours.  Considering  the  deployable  nature  of 
this  system,  a  fiber  coupled  scheme  would  reduce  or  eliminate  the  risks  associated  with 
misalignment  of  the  optical  board  as  well  as  the  time  required  to  change  the  alkali  being 
investigated. 

In  addition,  while  the  visible  alignment  laser  was  effectively  upgraded  prior  to 
this  experiment  to  enable  propagation  over  2  km,  efforts  at  longer  distances  would  likely 
benefit  from  a  new  configuration.  By  coupling  the  visible  alignment  laser  beam  to  the 
telescope  itself  such  that  the  beam  came  from  the  same  telescope  as  the  TDLAS  beam,  it 
would  be  possible  to  quickly  guarantee  alignment  over  longer  distances.  In  addition,  it 
would  decrease  the  time  required  to  coalign  the  visible  and  IR  beams,  a  task  that  must  be 
accomplished  at  every  new  distance  due  to  the  current  placement  of  the  visible  HeNe  on 
top  of  the  transmit  telescope. 

Several  issues  with  regards  to  platform  jitter  were  noted  at  the  2  km  range  of  this 
test.  With  a  baseline  noise  value  anywhere  from  0. 1  to  1  in  absorbance  at  2  km,  it  is  clear 
that  future  experiments  at  greater  path  lengths  will  require  a  significant  reduction  in  the 
effects  of  platform  jitter.  One  method  currently  being  investigated  is  the  use  of  a  1’  by  1’ 
retroreflecting  array  as  a  replacement  for  the  6  inch  mirror  currently  being  used.  First, 
the  increased  area  of  the  array  means  a  larger  effective  target  for  the  transmitted  beam, 
reducing  the  impact  of  telescope  jitter.  Also,  due  to  the  nature  of  a  retroreflector,  all  light 
incident  on  it  will  be  returned  along  its  initial  trajectory,  in  this  case  back  to  the 
transmitting  telescope.  This  property  would  significantly  decrease  the  influence  of  mirror 
misalignment  and  jitter. 
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In  terms  of  application,  given  the  necessary  upgrades,  similar  experiments  to  this 
one  could  be  conducted  at  even  greater  distances,  stretching  beyond  5  km.  Such  efforts 
would  be  valuable  in  continuing  to  evaluate  the  performance  of  alkali  laser  wavelengths 
in  the  atmosphere.  In  addition,  conducting  similar  tests  in  a  growing  number  of  locations 
would  enable  the  characterization  of  the  various  alkali  systems  in  the  environments  in 
which  they  would  be  used.  By  further  increasing  the  path  length  interrogated  and 
conducting  experiments  in  various  locations,  it  would  be  possible  to  more  accurately 
gauge  the  performance  of  each  alkali  system  in  a  given  high  energy  laser  scenario. 

Beyond  the  investigation  of  DPAL  wavelengths  in  the  atmosphere,  the  broad 
range  of  rotational  lines  interrogated  over  the  course  of  a  full  tunable  diode  laser  scan 
provides  an  opportunity  to  investigate  and  validate  the  broadening  rates  listed  in  the 
HITRAN  database.  While  this  possibility  was  investigated,  cleaner  data  would  likely 
provide  a  much  better  estimate  as  to  the  observed  broadening  rate  of  each  rotational  line. 
In  this  way,  the  rotational  dependence  of  the  broadening  rates  found  in  HITRAN  could  be 
investigated  with  greater  accuracy. 

Long  term  baseline  data  was  collected  over  the  course  of  this  experiment  which 
could  be  used  to  further  investigate  the  influence  of  turbulence  and  wind  speed  on 
baseline  fluctuations.  This  data  could  be  used  in  order  to  determine  whether  or  not 
turbulence  and  scintillometry  data  could  potentially  be  extracted  from  the  baseline  data  of 
the  spectra  collected  in  these  experiments. 
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Appendix  A 


Table  4:  HITRAN  08  data  for  the  lines  observed  in  the  890  to  910  nm  range 
surrounding  the  Cs  D1  emission  line. 
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Appendix  B 


Figure  13:  Atmospheric  absorption  spectrum  in  the  region 
surrounding  the  Cs  DPAL  emission  line,  taken  from  1545  to  1700  on 
17  October  2012.  The  scan  ranges  from  11,140  to  11,210  cm'1  and 
contains  ~30  absorption  line,  a  number  of  which  are  obscured  by 
baseline  noise. 
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Absorbance,  - (ln( It/Io)) 


i  1  i  1  i  1  i  1  i  1  i  1  i  1  i  1  i  1  r 


i  ■  i  i  i  ■  i  ■  i  .  i  ■  i  ■  i  ■  i  ■  i  ■  i 

12,480  12,485  12,490  12,495  12,500  12,505  12,510  12,515  12,520  12,525  12,530 


Wavenumber,  />,  (cm  1) 


Figure  14:  Atmospheric  absorption  spectra  from  12,480  to  12,530  cm  1 
collected  from  1628  to  1909  on  19  October  2012.  This  spectrum 
contains  ~25  weak  H2O  absorption  lines,  many  of  which  are  obscured 
by  the  noise  in  the  baseline. 
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Absorbance,  —  (ln(  It/Io)) 


Wavenumber,  (cm  l) 


Figure  15:  Atmospheric  absorption  spectra  from  12,554  to  12,570  cm  1 
collected  from  2159  on  19  October  2012  to  0050  on  20  October  2012. 
The  lower  noise  level  is  due  to  the  signal  averaging  used.  With  6  scans 
per  piezo  range,  the  signal  to  noise  ratio  was  improved  significantly. 
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Figure  16:  Atmospheric  absorption  spectra  in  the  Cs  DPAL  emission 
region,  from  11,150  to  11,210  cm'1,  collected  from  1905  to  2136  on  22 
October  2012.  This  spectrum  contains  ~30  strong  H2O  absorption 
lines,  some  of  which  are  obscured  by  the  noise  in  the  baseline,  as  well 
as  the  Cs  doublet  located  at  -11,182  cm'1. 
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Absorbance,  — 111(1/1$ ) 


12,555  12,560  12,565  12,570  12,575  12,580  12,585  12,590  12,595  12,600 

Wavenumber,  £%  (cm-1) 

Figure  17:  Sample  of  data  collected  in  the  Rb  region,  from  12,555  to 
12,600  cm'1,  on  23  October  2012  from  1515  to  1542  consisting  of  ~25 
weak  H2O  absorption  lines  as  well  as  the  Rb  Di  emission  lines  at 
-12,578cm1. 
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Absorbance,  -ln(I/Io) 


12,525  12,530  12,535  12,540  12,545  12,550 

Wavenumber,  t>,  (cm-1) 


Figure  18:  Atmospheric  absorption  spectra  from  12,525  to  12,550  cm  1 
collected  from  1551  to  1705  on  23  October  2012.  This  spectrum 
contains  over  30  weak  H2O  absorption  lines,  many  of  which  are 
obscured  by  the  noise  in  the  baseline. 
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